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Abstract 

By a systematic survey of the parameter space we confirm our surmise ^ that the 
Minimal Supersymmetric GUT(MSGUT) based on the 210 © 126 © 126 © 10 Higgs 
system is incompatible with the generic Type I and Type II seesaw mechanisms. The 
incompatibility of the Type II seesaw mechanism with this MSGUT is due to its 
generic extreme sub-dominance with respect to the Type I contribution. The Type I 
mechanism although dominant over Type II is itself unable to provide Neutrino masses 
larger than ~ 10“^ eV anywhere in the parameter space. Our Renormalization Group 
based analysis shows the origin of these difficulties to lie in a conflict between baryon 
stability and neutrino oscillation. The MSGUT completed with a 120-plet Higgs is 
the natural next to minimal candidate. We propose a scenario where the 120-plet 
collaborates with the 10-plet to fit the charged fermion masses. The freed 126-plet 
couplings can then give sub-dominant contributions to charged fermion masses and 
enhance the Type I seesaw masses sufficiently to provide a viable seesaw mechanism. 
We give formulae required to verify this scenario. 
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1 Introduction 


so(10) GUTs accomodate complete fermion families together with the superheavy right 
handed neutrinos (required by the seesaw mechanism[21 Ej) in a 16-plet spinorial irreps. 
They also provide very natural Higgs multiplets (126) (or 16 h x 16h) to generate the 
{U(1}b-l breaking) right handed neutrino Majorana masses and the {SU{2)l x U{1) 
breaking) triplet vevs required to implement the Type 1(2] and Type 110 mechanisms 
respectively. Thus models which incorporate high scale breaking of the SU{2)r x U{1)b-l 
components of the SO(IO) gauge symmetry provide an elegant and natural context to 
understand the clear indication of a seesaw connection between neutrino mass and GUT 
mass scales provided by the discovery of neutrino masses in the 50 milli-eV range (Am^ ~ 
2.2 X ~ ((10^GeU)^/(10^®GeU))^ ). In supersymmetric models with left right 

symmetric {SU{3) x SU{2)l x SU{2)b x U{1)b-l C G) gauge groups, R-parity {Rp = 
(^_)3(s-l)+ 2S^ ji] is a part of the Gauge symmetry and is naturally preserved IHl I71IH1 lUl ITH] 
till low energies in models using only B-L even vevs for the seesaw. They thus predict a 
stable LSP which is welcome as cosmological dark matter. With such a complement of 
virtues Supersymmetric SO(IO) GUTs of some variety or the other are considered lead¬ 
ing contenders by a majority of workers. The Supersymmetric SO (10) GUT based on the 
126 © 126 © 210 Higgs multiplets proposed long ago mma shares all these manifest 
virtues of LR supersymmetric models. It has the additional virtue of maximal simplicity 
from a parameter counting and representation economy point of viewjlll IT^ IT!?1 ITH] . It 
thus lays claim to the role of contributing the correct GUT symmetry breaking sector (AM 
Higgs) for a fully realistic Minimal Supersymmetric GUT (MSGUT). The other crucial 
component required to fully define a MSGUT is however the complement of fermion mass 
(FM)Higgs used to £t the charged fermion masses and mixings and the neutrino oscillation 
mass and mixing parameters which constitute the actual ‘data’ that a GUT must confront. 
SO(IO) permits an FM Higgs system based on 10 © 120 © 126 irreps. The first and third 
of these are have Yukawa couplings to the matter 16-plets which are symmetric under 
interchange of family indices and have been extensively considered. The 120-plet with its 
antisymmetric Yukawa couplings has so far been conceded much less importance, acting 
either as a perturbation to the main structure determined by the 10 and 126 plets miEii 
or else somewhat cursorily [nmniiiHiiin]. In this paper we shall develop a line of argument 
that culminates in a suggestion that a central role for the 120-plet is natural. 

In 1992, with LEP data in hand, Babu and Mohapatra^^l proposed that if only the 
10 © 126 irreps were used to fit the charged fermion masses and mixing then the matter 
fermion Yukawa couplings of the GUT would be completely determined. Hence the model 
would become predictive in the neutrino sector. The compatibility of this scenario with the 
(now better known) neutrino mass-mixing data has been the subject of an accelerating suc¬ 
cession of works since their proposal iiniiiniEilEaEiiEsi. These works have demonstrated 
the compatibility of the data with the generic (Georgi-Jarlskog [21] plus Type I and Type 
H seesaw) structural form of the fermion mass formulae dictated by SO(IO) clebsches and 
restricted FM Higgs content. Symmetric Yukawa couplings alone are quite successful but 
the possibility of improving the fits, particularly as regards the GKM GP phase, by using a 
subdominant contribution of the 120-plet has also been demonstrated |211III|. However the 
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generic fitting procedure does not fix the over all mass scale of the neutrino masses or the 
relative strength of the Type I and Type II masses. On the other hand these quantities are 
hxed in a fully specihed GUT model and it is thus crucial to investigate whether any given 
GUT model of this general type is compatible with both the generic hts and the actual 
neutrino mass data. In the 1990 ’s the construction of natural and fully consistent Minimal 
Left Right Supersymmetric models (MSLRMs) was accomplished [3 El IZI- Moreover these 
analyses showed clearly[2HlEI what was already noted at the beginningjTTlEZ] of the study 
of multi-scale Susy GUTs : that in LR Susy GUTs there are light multiplets which violate 
the conventional wisdom of the “survival principle”. Thus it was clear P EHj that Susy 
GUTs required the use of calculated rather than merely (“survival principle”) estimated 
masses for RG analyses. The RG analysis carried out in this work already indicated that 
the use of calculated spectra forces together the various possible intermediate scales into 
a narrow range close to the GUT scale resulting in an effective “SU(5) conspiracy” i.e the 
necessity for single step breaking of SO(IO). 

A calculation of the full GUT spectrum and couplings in various SO(IO) GUT models 
required complete knowledge of the ’’Glebsches” relating S'O(IO) and Gsm x U{1)r group 
labels on SO(IO) helds - at least for tensorial irreps - before progress in the Renormalization 
Group (RG) analysis could be made. Furthermore the compatibility of fermion hts with 
GUTs could also be investigated only if Glebsches were also available for the spinorial 
(16-plet) irrep. The requisite computations (based on decomposition of SO(IO) labels into 
those of the SU{4:) x SU{2)l x SU{2)r (“Pati-Salam”) sub-group ) were presented by us 
in ra. Using these Glebsches hrst the mass matrix of the MSSM type doublets ([1, 2, ±1]) 
and classic GUT AR ^ 0 process mediating triplets ( [3,1, ±|] [201 and then a complete 
calculation[2I], of all the couplings and mass matrices of the MSGUT was given by us. 
With similar motivations two calculations, one in parallel and cross checking with oursj32j. 
and another [221 quite separate, which both used the same somewhat abstract |22] method 
(but different phase conventions) to calculate the “Glebsches” of tensor (but not spinor) 
multiplets, appeared. In j2ni cm we provided the complete (chiral and gauge) spectra, 
neutrino mass matrices, gauge and chiral couplings and the effective d = 5 operators for 
Baryon violation in terms of GUT parameters and MSSM helds. Thus the stage was 
laid both for a completely explicit RG based analysis of the MSGUT and a completely 
specihed investigation of the compatibility of Type I and Type II seesaw mechanisms. The 
calculation of threshold ehects based on these spectra was described inj2Il C] and shown 
to controvert the conventional wisdom [22] which doubted the stability of the successful 
unihcation (and mt,sin^6'vi/ prediction |221) of the one-loop MSSM coupling evolution [HHl 

I2I1I2HI. 

The question of the htting of the charged fermion masses, the relative size of the neutrino 
mass splittings and the quark and lepton mixing matrices had been extensively analyzed |13l 
linillllimimEllESI using the generic formulae valid in SO(IO) GUTs with only|T2| 10,126 
FM Higgs representations. However this procedure does not constrain the relative size of 
the Type I and Type H contributions or the overall scale of the neutrino masses. Within 
a fully specihed theory like the MSGUT, however, both these parameters are specihed. In 
|22| the question of assuring Type H (over Type I) dominance on the basis of the Yukawa 
values found in the successful generic Type H hts |21[ |22j was considered. The authors 
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concluded on the basis of order of magnitude estimates that the Type I mass would emerge 
as too large unless the B — L breaking scale was raised. Furthermore they found that the 
mass of the SU{2)l triplet held whose tadpole vev gives rise to the Type II seesaw must 
be lowered to enhance Type II seesaw masses. To ensure these they considered extending 
the model by modifying the GUT scale breaking to be via S'f/(5)(this required introducing 
a 54 — plet ) or else by perturbing it by introducing a 120-plet. 

Our point of view, program and results on these questions are related but still quite 
distinct from the estimates and scenarios of 123- In[I] we took up the question of the 
compatibility of the MSGUT with the successful generic 10 + 126 FM Higgs hts. We 
surveyed the variation of the relative strength of Type I and Type II seesaw mechanisms, 
as well as their absolute magnitudes while keeping in view the viability of the relevant 
regions of the GUT parameter space i.e with the stability of the basic unihcation scenario. 
Our graphical survey yielded the striking observations that : 

• For generic as well as special values of MSGUT couplings and FM Yukawas taken 
from the generic Type II seesaw fit, the Type II seeaw itself implies that it is highly 
subdominant to Type I seesaw in the MSGUT unless the relative strength of Type I 
versus Type II is adjusted to a very small value. The primary reason for this relative 
dominance of Type I is not the value of the symmetry breaking scale but rather the 
value of the effective coupling in the Type I seesaw formula implied by the Yukawas 
found in the Type II generic fit (which simply assumes Type II is dominant) 

• The maximal values of the Type I seesaw masses attainable in the MSGUT are at least 
one order of magnitude short of those reguired by atmospheric neutrino oscillations . 

Note that on both accounts our conclusions are at variance with the j2S| although they 
also hnd difficulties in ensuring Type II dominance over Type I. 

In the present paper we present a detailed survey of the parameter space to dehnitely 
conhrm the above conclusions and thus pose a stringent challenge to the viability of the 
original proposal of Babu and Mohapatra^H] in the context of the MSGUT. Already in 
we noted that the variation of the unihcation stability monitoring parameters (USMPs)(i.e 
A{LogMx), I^{Sin^9w), with the fast control parameter ^ of the MSGUT exhib¬ 

ited very striking sharp peaks and dips (with the characteristic appearance of poles of a 
function of f). In many cases these spikes were at the values of points of increased symme¬ 
try uncovered by the analytic solution of the AM spontaneous symmetry breaking problem 
in the MSGUT (TUI ED E2j- addition there were other spikes (not obviously related to 
points of extended symmetry) whose investigation was called for by the “tomogram” of the 
MSGUT parameter space that we had provided ED. InjUUEDED an analytic parametriza- 
tion of the SSB in the MSGUT in terms of the solutions of a cubic equation (linear in 
a fast control parameter was achieved. In particular in |1()1 ED ED fhe parameter ^ 
was eliminated -using the cubic equation- for the solution values x and thus a very di¬ 
rect,elegant and unihed parametrization of the pole and zero structures characterizing the 
MSGUT AM-SSB was achieved in terms of the variation of a single complex variable x. 
This allowed the identihcation of additional special points. In particular, motivated by our 
observations [D, of a set of points associated with possible growth of the Type I and Type 
II seesaw contributions [101 ED ED- Such growth had already been previously observed at 
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some points in P] but found to give unviable USMPs. Thus to complete the program otni 
we have now investigated the behaviour of the Type I and Type II coefficient functions 
dehned by the precise clebsches and fermion mass matrices of the MSGUTinUID over the 
complex X ( equivalently to ^ ) plane : including all the special points described above. The 
loopholes offeredPlIlSj by the special points are illusory rather than real : typically some or 
all of the USMPs explode due to the very growth which is being invoked to strengthen one 
or the other seesaw. Our conclusions are thus unchanged from those of m the MSGUT 
has suffered a failure of its most critical features and is thus now defunct. 

A way out of this impasse appears by allowing in the 120-plet which was earlier dis¬ 
criminated against on grounds of convenience and simplicity alone. ‘Lamppost logic’ may 
have failed here as so often before. Nature does not heed the wishful thinking behind such 
logic: the morning often reveals that the keys sought for all night lay in the gutter just 
beyond the lamppost ! Another alternative is that SSB in the MSGUT may actually be 
controlled isniiini by the explosive growth of the gauge coupling above Mx due to gaugino 
condensation in the coset S'O(10)/Gi23 which drivesjin] chiral condensates of the AM Higgs 
helds whose value may be calculable due to the power of supersymmetric holomorphy. In 
such a scenario the crucial composition of the zero mode Higgs doublets must be evaluated 
without using the “cubic” solution mentioned above. 

In Section 2, referring the reader to the original papers Ennsnin for all derivations, 
we present the basic formulae jIHL governing the Type I and H seesaw mechanisms in 
the MSGUT and identify the Seesaw Monitoring Parameters (SMPs) {R, Fj, Fu) that we 
will use for our graphical presentation of the seesaw relevant topography of the parameter 
space. We list the obvious points of enhanced gauge symmetry. This is accompanied by 
an Appendix where these formulae are also evaluated in the (Rational) function (of x) 
form advocated by|lH]. This allows the identihcation of an additional list of exceptional 
points requiring special treatment. In Section 3 we describe the analytic and graphical 
investigation of the behaviour of the USMPs and SMPs as we scan over the complex x 
plane for generic and for exceptionally favourable values of the ‘slow’ AM Higgs param¬ 
eters EH i.e, X (equivalently . By covering the behaviour at the enhanced symmetry 
points and exceptional seesaw points (collectively called ESPs)and scanning the whole x- 
plane we eliminate all escape routes and pin down the MSGUT to its doom. In Section 
4 we propose that the inclusion of the 120-plet in a very specihc role can solve the prob¬ 
lem with the small neutrino masses found for the MSGUT and give the fermion mass 
formulae for the NMSGUT along with the doublet Higgs mass matrix Ti. This specihes 
the composition of the massless doublets once detTi = 0 is imposed. Thus the stage is 
set for the accomplishment of a program analogous to that accomplished for the MSGUT 
[iiiiniisniiiniEiiEaEsiiiiiiaiiiiin!, which can be implemented without ambiguity. We 
conclude with a discussion of the outlook for future work. 
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2 Seesaw formulae and SMPs for the MSGUT 


2.1 MSGUT couplings, vevs and masses 

The MSGUT is the renormalizable globally supersymmetric S'O(IO) GUT whose Higgs 
chiral supermultiplets consist of Adjoint Multiplet (AM) type totally antisymmetric ten¬ 
sors : 210($ijfcz), 126(SijfcZm), 126 = 1...10) which break the GUT symmetry 
to the MSSM, together with Fermion mass (FM) Higgs lO-plet(Hj). The 126 plays a 
dual or AM-FM role since it also enables the generation of realistic charged fermion and 
neutrino masses and mixings (via the Type I and/or Type H mechanisms); three 16-plets 
^a(A = 1,2,3) contain the matter including the three conjugate neutrinos (h^). 

The superpotential (see [121 EDI EH E2| for comprehensive details ) contains the mass 
parameters 


m : 210^ ; M : 126 ■ 126; Mh : 10^ (1) 

and trilinear couplings 

A : 210=^ ; 7] : 210 ■ 126 TM; 7 © 7 ; 10 • 210 ■ (126 © TM) (2) 

The GUT scale vevs that break the gauge symmetry down to the SM symmetry (in the 
notation ofED]) are fTTlfT^ 


{(15.1.1))2io 

(0afecd) 2 ^ctbcdef^ef 

(3) 

{(15, 1.3)>2i„ 

i^aba/s) ^ 

(4) 

((1, I5 1))210 

i^apps) = P^apps 

(5) 

((10,1,3))t^ 

(^iaSso) ~ ^ 

(6) 

((I0,1,3))i26 

(^ 24679 ) ~ 

(7) 


The vanishing of the D-terms of the SO(IO) gauge sector potential imposes only the con¬ 
dition |(t| = |a|. Except for the simpler cases corresponding to enhanced unbroken gauge 
symmetry (S'U(5) x U{1), SU{5), 03 ^ 2 , 2 ,b-l, G 3 , 2 ,r,b-l etc)[TniE21, this system of equations 
is essentially cubic and can be reduced to the single equation PI for a variable x = — Acn/m, 
in terms of which the vevs a,uj,p, a, a are specihed : 

8x^ - 15x^ + 14a; - 3 = -^(1 - xf (8) 

where ^ Then the dimensionless vevs in units of (m/A) are uj = —x ^01 and 

(a;^ + 2a; — 1) _ a;(5a;^ — 1) 2 Aa;(l — 3a;)(l + a;^) , . 

“= (l-x) • ( 1 -xP • = - 0^75- 

This exhibits the crucial importance of the parameters x. Note that one can trade pm 
E21I1S] the parameter ^ for x with advantage (using equation (jH))) since ^ is uniquely hxed 
given X. By a survey of the behaviour of the theory as a function of the complex parameter 
X we are simultaneously covering the behaviour of the three different solutions possible for 
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each complex value of We thus change our graphical presentation in terms of three plots 
versus einma to a single plot versus x for each USMP or SMP. Moreover, as emphasized 
bv|4.fj analysis of the spikes observed in (sum is also facilitated by the use of the parameter 
X. Using the above vevs and the methods of HDl we calculated the complete gauge and 
chiral multiplet GUT scale spectra and couplings for the 52 different MSSM multiplet sets 
falling into 26 different MSSM multiplet types (prompting a natural alphabetization of 
their naming convention [30j !) of which 18 are unmixed while the other 8 types occur in 
multiple copies. The (full details of these) spectra may be found inj^HlEIl and equivalent 
results(with slightly differing conventions) are presented injHSj. A related calculation with 
very different conventions has been reported in [SSI mi. The initially controversial relation 
between the overlapping parts of these papers was discussed and resolved in jH] . 

Among the mass matrices is the all important 4x4 Higgs doublet mass matrix jsniEii 

n : 


( -Mh 

-7\/3(i.j + a) 

7\/3(cj — a ) 

V - 0-7 


+7\/3(ii) - a) 

0 

— (2M + 4r}(a — to)) 
— 2r}a\/3 


-'yVsiiU + a) 

— (2M + 4r)(a + u))) 
0 
0 


—'ja 

0 

— 27]a\/3 

—2m + 6A(a4 — a) 


( 10 ) 


Ti can be diagonalized by a bi-unitary transformation |TQ1IS21E|: from the 4 pairs of Higgs 
doublets arising from the SO(IO) helds to a new set of helds in terms of 

which the doublet mass terms are diagonal. 

U^TiU = ....) 

( 11 ) 


To keep one pair of these doublets light one tunes Mh so that DetTi = 0. This matrix 
can then be diagonalized by a bi-unitary transformation yielding thereby the coefficients 
describing the proportion of the doublet helds in the 10,126,126, 210 GUT multiplets 
present in the light doublets : which proportions are important for many phenomena. 
In the effective theory at low energies the GUT Higgs doublets are present in the 

massless doublets in a proportion determined by the hrst columns of the matrices 

U,U : 

E«Mx : ; ai = Uh 

^ a.ij(i) ; ai = Ua (12) 

The all important normalized 4-tuples a, d can be easily determined by solving the zero 
mode conditions: Ha = 0 ; dP^H = 0. 


2.2 RG Analysis and USMPs 

In jSH C] we discussed at length the use of plots of the USMPs versus the fast parameter 
^ to investigate the question : 

Are the one loop values of adMx), Sin‘^9w and Mx generieally stable against superheavy 
threshold corrections ?. 

We followed the approach of Hall|13] in which the mass (Mx = \JA\d + + 2\p + ojP) 
of the baryon number violating superheavy ([3,2,±|] or X-type) gauge bosons is chosen 


6 





as the transition scale between the effective MSSM and the fnll SO(IO) GUT with all 
snperheavy fields retained. Thns Mx is nsed as common physical snperheavy matching 
point {Mi = My = Mx) in the eqnations relating the MSSM conplings to the SO(IO) 
conpling : 


1 

cti{Ms) 


M 


+ 87ibiln"-^ + - 4:7iXi{Mx) 


aG{M_ 


X 


M: 


(13) 


See jlSl for details. In this approach it is recognized that above the scale Mx the 
effective theory changes from the MSSM to Susy S'(9(10) model structured by the com¬ 
plex superheavy spectra which we have computed so that the theory appears as unbroken 
SO(IO) only to leading order. Thus we compute the corrections to the three USMPs 
LogioMx, sin‘^9w{Ms), aQ^{Mx) as a function of the MSGUT parameters and the answer 
to the question of stability of the perturbative unification is determined by the ranges of 
GUT parameters where these corrections are consistent with the known or surmised data 
on LoqiqMx, sin‘^9w{Ms), ac- The consistency requirement that the SO(IO) theory remain 
perturbative after threshold and two loop correction and ac not decrease so much as to 
invalidate the neglect of one-loop effects in the couplings A, rj, 7,7 is also appropriate. We 
implement these by the requirements that the USMPs obey 


AgI 

= |A(a7(AG))l < 10 


Ax 

= A{LogioMx) > — 1 


Aw 

= \A{sin'^9w{Ms))\ < .02 

(14) 


See jSH Cl for a detailed discussion. We find the corrections 

M' 

A^^^\LogioMx) = .0217 + .0167(56; + 36; - %h'^)Logw^ 

Mx 

M' 

A^^'^\sin^9w{Ms)) = .00004 - .00024(46; - 9.66; + bM'^)LogiQ^ 

Mx 

M' 

A^^^\aQ^{Mx)) = .1565 +.01832(56;+ 36; + 126'3)Loc/io7^ (15) 

Mx 

Where 6 ' = 167 r^ 6 ' are 1-loop fd function coefficients (/9j = bigf) for multiplets with mass 
M' (a sum over representations is implicit). Note a minor, but crucial correction, in the 
above equations relative to EUd: there and here, we actually used as the argument 
of the logarithms in these formulae. This is the value given by the algebra, which is then 
conventionally approximated by since the difference is of order gm- However in our 
case we can retain the “exact” form and still calculate the threshold corrections since those 
depend on ratios of masses. This improves the analysis and makes it more convenient. 
Moreover (see below) it means that the overall scale parameter that we use , namely the 
coefficient m of 210^ in the superpotential, is not fixed at the one loop unification scale 
10i6.25(7g^ but slides with the corrections to Mx- These corrections are to be added to 
the one loop values corresponding to the successful gauge unification of the MSSM : Using 
the values 


a%{Mx)~^ = 25.6 ; GeV ; Ms = 1 TeV 

a^\Ms) = 57.45 ; a^^Ms) = 30.8 ; a^^{Ms) = 11.04 (16) 
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-.08 


A‘^-^°°P{sin‘^ewiMs)) = .0026 


the two loop corrections are 

/^2-ioop^-i^Mx) = -.546 (17) 


We see that in comparison with the large threshold effects that might be expected 185] in 
view of the large nnmber (504) of heavy fields the 2 loop corrections are qnite .small 1811. 
The striking new resnlt of EH was thns the explicit demonstration that the coniectnrejd5j 
concerning the “Fntility of high precision SO(IO) calcnlations” was in fact false and that 
a proper calculation based on calculated spectra yields well defined and consistent results 
for significant regions of the MSGUT parameter space (see below). 

The parameter ^ = XM/rjirt is the only numerical parameter that enters into the cubic 
eqn.(jH)) that determines the parameter x in terms of which all the superheavy vevs are 
given. It is thus the most crucial determinant of the mass spectrum . The dependence 
of the threshold corrections on the parameters A, rj, 7,7 is comparatively mild except when 
coherent e.g when many masses are lowered together leading to ac explosion, LogMx 
collapse or large changes in sin^ 0w{Ms). 

The typical behaviours of the USMPs can be seen in Figs 1-6. 

The parameter ratio m/X can be extracted as the overall scale of the vevs. Since the 
threshold corrections we calculate are dependent only on (logarithms of) ratios of masses, 
the parameter m is fixed in terms of the mass My = Mx of the lightest superheavy vector 
particles mediating proton decay by : 


Ax 

\m\ 


A(io9.o^) 

^q16.25+Ax GeV 

g\J^\a + w\‘^ + 2\p + (jj\^ 


( 18 ) 


The presence of the factor 10'^^ in the formulae for the neutrino masses will play a crucial 
role as one of the hands of the “scissor” operated by Baryon Decay and constraints and 
neutrino oscillation data. The lowering of the “slow diagonal” parameters A, rj tends to make 
fields light and thus give large negative corrections to Ac (as can be seen in the expansion of 
the darker areas in Fig.4 relative to Fig.3). Such growth of the gauge coupling can invalidate 
the self consistency of the perturbative approximation used through out. Conversely when 
ac decreases too much the neglect of Yukawa loops becomes moot. Thus A, g are effectively 
restricted to magnitudes of order 1. Note, however, that A^y, Ax are insensitive to the 
precise values - whether real or complex - of A, 7 as long as their magnitudes are ~ 1 : as 
can be seen in the pairs Fig.l and Fig. 6 and Fig. 2 and Fig. 5. In addition there are the 
“slow off diagional” parameters 7,7 whose effect is quite mild, and can be taken to be any 
values ~.l to 1 ., without any appreciable change in the numerical plots. 


2.3 Fermion mass formulae and SMPs 

The Dirac masses of matter fermions in the MSGUT at scale Mx are 

= v{h + f) 
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lm[X] 31 ^ e^^p|(jt(Cont->-.0!2-.01 ,X>^,.02),a= 1.0, 17 = 0.95, j= 0.7,gambar=0.6 > 



-3 -1.5 0 1.5 3 


Figure 1: Contour Plot of the threshold corrections to Siri^Oyj on the complex x plane. 
Contours at /S.w = —-02, —.01, .01, .02. The shading progresses from black (< —.02) to 
white > .02. 


m!' = v{h — 3/) 

= v{rih + r2f) 

m} = v{rih — 3r2f) (19) 


Here we have extracted v = 174:GeV. Note the characteristic (Georgi-Jarlskog) j2E] factors 
of (-3) in the lepton masses relative to the quark masses. Here [201 Ell Cl 


h 


ri 


2v^hQ!i sin/?; / = —4y-i/a 2 sin/? 

— cot /?; r 2 = — cot /? 

(y.\ 0-2 


( 20 ) 


and (h, /) are the couplings in the MSGUT superpotential. They are 3x3 symmetric 
matrices with rows and columns indexed by the SO(IO) family indices A,B = 1,2,3 |101 
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lm[X] 


i){Plot(Cont->-1-.5,.5,1>,a= 1.0,1?= 0J95, j= 0.7 flam bar= 0.6 > 



R*[X] 


Figure 2: Contour Plot of the threshold corrections Ax to Log\i^{Mxl^GeV^ on the 
complex X plane. Contours at Ax = —1, —.5, .5,1. The shading progresses from black 
(< —1) to white > 1. 


EDI EB HI- Notice that since |aj| < 1 by unitarity, and \h,f\ < .5 by perturbativity, 
the inclusion of ai,a 2 in the definitions can yield subtleties only at zeros of ai,a 2 - The 
Majorana mass parameters My,Mp of the left and right handed neutrinos (defined as the 
coefficients of z/i/, h h in the superpotential are^Tl IT] 

Mp = fa 

M, = rj ( 21 ) 


where jnil IB 


. ia\/3 

02 sin (3 

Mq = 2(M + p(3a —p)) 


( 22 ) 
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lm[X] 


|islPlot(Cont->-10-5,5,10>,(:A= 1.0,7= 0.95, y= 0.7,gambar=0.6 > 



^----- Re[X] 

-3 -1.5 0 1.5 3 

Figure 3: Contour Plot of the threshold corrections Ac to Uq^^Mx) on the complex x 

plane. Contours at = -10,-5,5,10. The shading progresses from black (< —10) to 
white > 10 . 


ra = -2i'/3{ai'y+ 2y/3r]a2){—){^) sin f3 

02 Mo 

From these one obtainsjSI] the Type I[2]and Type II seesawj^ Majorana masses (in the 
conventional normalization W = M^vvI2 + ...)of the low energy neutrinos: 

Ml = vr^h ; = 2vrj,f 

ia 2 sin /3 v 


Ti = 


2v/3 d 
n = (h - 3/)/“^(h - 3/) 


(23) 


For investigation of the Type I seesaw we will find it useful to define functions F,,Fn by 
((1.70 X 10 ^eV) = v'^/Mx and we have omitted irrelevant phases) 

Ml = (1.70 X 10“^el/) sin/3F7 h 
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lm[X] 


|iGlPloKCont->-10,-5,5,10>A=0.1, 7= 0.7, y= 0jS,gambar=0.S > 



Figure 4: Contour Plot of the threshold corrections to aQ^{Mx) on the complex x plane 
illustrating the drastic effect of decreasing the value of the diagonal parameter A. Contours 
at Ag = —10, —5, 5,10 the shading progresses from black (< —10) to white > 10. 



12 



lm[X] 


i){PlotCCont->-1,-.5,.5,1>j(A=0.1, ri= 0.7, y= 0.6,gambar=0B > 


3 



^----- Re[X] 

-3 -1.5 0 1.5 3 

Figure 5: Contour Plot of the threshold corrections Ax to Logio{Mx/lGeV) on the 
complex X plane showing the minor effect of reducing the value of the diagonal parameter 
A. Contours at Ax = —1,—.5,.5,1 the shading progresses from black (< —1) to white. 
> 1 


In these formulae the scale parameter m has been eliminated in favour of Ax using eqn (j18jl . 
Explicit forms for Fj, Fu in terms of the parameter x are given in the Appendix. We 
also define the ratio R = \Fi/Fn\ to capture the relative strength of the two seesaw 
contributions. A crucial point is that if one wishes to maintain perturbativity and work 
with GUT Yukawa couplings < 1 then the large tan/? scenario is necessary. The Type 
I and Type II fits of the fermion mass-mixing data are carried out 1213 El 122 El EH] 
assuming only the characteristic form of eqns. mm to derive “sum rules” [HEniElEI] 
among the different mass matrices at the GUT scale. In these generic fits, the freedom to 
choose the coefficients ri, r 2 , ra, r 4 is assumed to be compatible with the theory in which the 
mechanisms are embedded. The hts then yield the neutrino masses only up to an overall 
scale and the required relative strength of Type I and Type II is simply assumed. In pQ 
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lm[X] 31 ^ e^^PloKCont->-j02,-.01,.01 ,j02>j(A=0.1, r}= 0.7, y= 0.6, gambar=0.S > 



-3 -1.5 0 1.5 3 

Figure 6: Contour Plot of the threshold corrections to Siri^Ow on the complex x plane 
showing the minor effect of reducing the value of the diagonal parameter A. Contours at 
Aiy = —.02, —.01, .01, .02 the shading progresses from black (< —.02) to white > .02. 


we showed that these assumptions are unjustified for the MSGUT which seems to impose 
extreme dominance of Type I over Type II hts and moreover is incapable of achieving 
neutrino masses as large as .OSeGas required by experiment. Our argument was based on 
the following observations : 

• When a pure Type II fit is assumed one finds that 1221 El ESI the maximal value of / 
eigenvalues is ~ 10“^ while the corresponding values for h are about 10^ times larger. 
This is a structural feature due to the approximate equality ofmb{Mx) and mr{Mx) 
which ensures that no Georgi-Jarlskog type contribution can he important for the third 
generation. As a result n ~ 10^/. This implies that the ratio R dehned above must 
obey R < 10“^ for the pure Type II not to be overwhelmed by the Type I values it 
implies. Our initial survey of the MSGUT parameter space showed that, generically, 
this did not happen and pure Type II could not work. Even near exceptional points 
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no escape was allowed by the need to maintain viable USMPs. In the next section we 
describe a much more exhaustive survey of the MSGUT to support this conclusion. 
Note that this is the crux of the reason for the failure of the Type II seesaw and it 
applies equally to “mixed Type I- Type II ” hts^Hj. 

• The Type I hts [THl OH] yield, typically, n 5/ 3 and thus we see that the 

magnitude of the function Fj will need to be greater than about 10^ in order that 
neutrino masses as large as the .05 eV required for compatibility with neutrino mass 
squared splittings as large as A(m '^)23 ~ .002eld^ be achievable. This was seen to be 
generically un-achievable inpp. In the next section we prove that this requirement 
combined with that of viable USMPs makes it impossible for the Type I £t to reach 
viability. This will be interpreted as indicating a missing piece to the puzzle, the 
omission of the 120-plet. 

The estimate of the magnitude of n chosen is obviously critical to the disqualihcation 
of the MSGUT Type I masses. A value value h ~ 10 would drastically modify our 
conclusions. The 3 generation fits performed so far are all simply numerical -although 
guided by analytic insight into where in the parameter space a £t might occur [13 ESI- 
Therefore every £t will have its own value of the eigenvalues of nnf which might in 
principle have widely different magnitudes. However no such exceptional Type I £t 
has been found so far. Although an analytical proof of this is still elusive - because 
of the very complexity that obstructs the determination of an analytic £t, there 
are heuristic arguments why a large h £t is unlikely to exist. The charged fermion 
£t succeeds because the second generation couplings of the 126 dominate those of 
the 10 so that the Georgi- Jarskog mechanism can operate. Although the opposite 
(^33 >> / 33 ) is true for the third generation so that no such mechanism is effective 
(as it should not since at Mx the b, r masses are roughly equal) yet the the values of 
/ in the 2 — 3 sector are not so small that the eigenvalues of h become larger than 1. 
The Type I neutrino masses are found to obey a normal hierarchy where the second 
and third generation neutrino masses are much larger than the first generation so 
that the 2 — 3 block of n is dominant over the rest of n. The large neutrino mixing in 
the 2 — 3 sector implies that the matrix elements n in that block are roughly of the 
same order of magnitude. These elements are generically found to be of magnitude 
less than 1 : the inverse hierarchical structure of / ensures that the normal hierarchy 
in h — 3/ does not get amplihed by the f~^ factor in the Type I seesaw formula to 
the extent of increasing the maximal eigenvalue of h beyond .5 or so. An analytic 
proof of this is still lacking. Thus in principle there is a loophole which is still 
unclosed : exceptional Type I fits with maximal eigenvalue greater than 10 , if found, 
will need to be reanalyzed separately. We emphasize that the question is not one of 
hnding improved or optimal hts but only of the order of magnitude of the maximal 
eigenvalue of h. We shall assume that such exceptional fits do not exist. 
In other words we assume that the Type I fits of j 151 [25^ are generic 
and representative as regards the order of magnitude of h. It is important 
to build a data base or other more compact description of all passable achieved hts 
that may be used to survey all the possibilities. Thus a naive estimate for the typical 
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^ hp- 

maximal magnitude of h : Max n ~ ^ >> 1 is actually far off the mark. As 
explained above the characteristic off diagonal structure of the successful Type I hts 
in the generic Babu-Mohapatra proposal, seems to be intimately bound up with this 
difference between the naive and actual magnitudes of (the eigenvalues of) fi. All 
these remarks are illustrated by the values taken from the published example in ra ( 
which is, however, not a completely satisfactory £t since the electron mass is fit only 
to an accuracy of about 10%). 


h 


j -.00016 - .00013* 
-.00046 - .00035* 
\ -.00642 + .002i 


-.000456 - .00035* 
-.00584-0.00101 
.03142 


-.00642 + .002i 
0.03142 

-.55633 + .095531* 


(26) 


/ 


/ 0.00023 + 0.00018* 
0.00018 + 0.00017* 
V 0.00285 - 0.00053* 


0.00018 + 0.00017* 
0.00697+ 0.00164i 
-0.01347-0.00158* 


0.00285 -0.00053* \ 
-0.01347- 0.00158* 
0.00053- 0.037858* ) 


(27) 



/ 2997.76+ 347.393* 
115.886 - 418.543* 
V 125.607- 173.883* 


115.886-418.543* 
51.4218 + 16.1651* 
-40.8494 + 15.3566* 


125.607- 173.883* \ 
-40.8494 + 15.3566* 
-18.5738 + 4.48526 *) 


(28) 


ntz = Chs - 3 / 33 ) 7/33 = 6.37485 + 6.98329* (29) 


/ .00278 + .00284i .00712 + .00287* 

.00712 + .00287* .09225 + .02520* 

V.02697+.00393* -.12975 -0.08165* 


.02697+.00393* \ 
-.12975 - 0.08165* 
0.00256 + 0.15193* / 


(Eigenvalues [h.h7 2 = {0.2771,0.06298,0.01088} 


(30) 


(31) 


The magnitude of the maximum eigenvalue of n is seen to be almost 35 times smaller 
than the “generic argument” one! Thus generic arguments for this magnitude need 
to be taken with a large lump of salt !. 


3 USMP and SMP plots : Gloomy Pictures ? 

We wish to survey the USMP and SMP variation over the complex x plane. We will use 
generic but hxed values with magnitudes ~ 1 for the ‘diagonal’ slow parameters [HU C] 
(A,*]}. We have checked that for any such (real or complex) values there is no signihcant 
modihcation of our conclusions. If one lowers A, rj below about .2 or so the USMPs exhibit 
pathologies due to mass scale lowering. Specihcally, an increase of Ac beyond a value 
of 10 or a decrease below (say) —10 makes the perturbative formulae used throughout 
unreliable. Similarly varying the ‘non-diagonal’ slow parameters 7,7 to small values, though 
not forbidden by USMPs IT] is also infructuous for purposes of viable SMPs. 
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X 


Symmetry / remarks 

1/2 

-5 

SU{5) ,zero of p2 

-1 

10 

SU{5) flipped,zero of p 2 

0 

3 

SU{3) X SU{2)l X SU{2)R x U{1)b-l 

1/3 

-2/3 

Flipped SU{5) X U{1) 

±i 

3^ 6 i 

SU{3) X SU{2)l X f/(l)R X U{1)b-l 

0.123765 

1.93008 

zero of ps 

0.646451 - 0.505389 i 

4.86831 + 3.29252 i 

zero of p 3 

0.646451 + 0.505389 i 

4.86831 - 3.29252 i 

zero of p 3 

{3±iV7)/8 


zero of q 2 

-3.46301 

28.2958 

zero of ps 

0.262212 ± 0.388123 i 

2.13469 T 3.38025 i 

zero of ps 

0.358184 ± 0.133971 i 

-0.449266 T 2.30098 i 

zero of ps 

0.1984 

1.167 

zero of ^3 

-0.099 ± 2.24 i 

1.596 T 15.572 

zero of gg 


Table 1: Special points with extended gauge symmetry and/or candidate seesaw enhance¬ 
ment. 

Besides the survey across generic regions of the complex x ( equivalently the 3 ^ ) plane(s) we 
must also examine the possibility that the behaviour at certain exceptional 

points may invalidate the generic trend of seesaw exclusion. For this purpose we provide in 
Table I a list of the exceptional points . These consist on the one hand of a list of points of 
enhanced gauge symmetry where the USMPs typically show spikesjSIj and on the other of 
special points where the SMPs may exhibit spikesjUHB]. The complete analytic pi U HH] 
identihcation of the latter set of points using the convenient parametrization of vevs in 
terms of Rational functions of the variable x alone (i.e with ^ eliminated in favour of x) is 
related in the Appendix. 

3.1 The fatal weakness of Type II 

Let us begin with a discussion of the viability of the Type II mechanism. The SMP R we 
proposed PP to monitor the strength of Type I relative to the Type II seesaw is particularly 
simple to analyze. It is the magnitude of a rational analytic function ( the ratio of a 
polynomial of degree 10 to a polynomial of degree 9) of x and is dependent on no other 
parameters. 


Fi Ax - l){/x - l)q 2 q'i .oon 

^Fj/ U{Ax-l){x^+ l)qi' 'Q{x)^' ^ ^ 

(the polynomials g'are dehned in Table II in the appendix) 

Therefore R is bounded from below by 0 and its minimae are isolated : they are the zeros 
of P{x). As |a;| —>■ cxo it grows as |a;| (with coefficient 3/64). For |a;| > 3 ^ i? > 10“^ or 
so. This is clearly visible in the 3D plot shown in Fig. 7. Thus only for |x| < 3 do we 
need to even consider the possibility of dominant Type II seesaw. This is in fact the region 
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which contains the zeros of R namely x = {1/3,1, {(3 ± z-\/7)/8}, {0.198437, —.0992186 ± 
2.24266*}} of which the latter two sets are the zeros of q 2 and gg respectively. A picture of 
the variation of R in this region of the x plane is given by the 3D plots of R~^ in Fig. 8 
which clearly shows that, as expected, the zeros of the numerator dominate the behaviour. 



Figure 7: 3D Plot of the relative strength R of Type I and Type II seesaw mechanisms on 
the complex x plane showing that the region |x| > 3 is irrelevant for promoting Type II 
dominance. 


The real zeros of R show up as “cliffs” (instead of narrow peaks ) due to periodicity in 
Arg X which slices them vertically, while “ridges” observed in these graphs are associated 
with the complex zeros of ^ 2 ,^ 3 ) (one needs to magnify the graph to see the ridge form of 
the zeros of ^ 2 )- There is a strong preference for a narrow range of Arg x for minimizing R 
(as also seen directly in the plots of R in the neighbourhood of these zeros discussed below). 
Thus the growth of R~^ is seen to be well localized to the regions around the zeros of R. 
From this survey as well as the theory of analytic functions it is clear that the behaviour 
of R and the USMPs in the neighbourhood of its zeros will be sufficient to decide whether 
any Type II solutions are viable. We now discuss this in greater detail. 

Since Type II dominance desires as small a value as possible of R we need to approach 
as closely as viable to the zeros of P{x). Among the zeros of R, however, are also the 
zeros of 02 : {x = {1/3,1}, {0.1984,—.099 ± 2.24*} (zeros of The decrease of \a 2 \ 
is limited by the requirement of maintaining perturbativity in the superpotential Yukawa 
coupling ~ 3//a2 < 1. Type II dominant fits yield/ about 10“^ j2Zl EH 123 C] • Even if 
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Figure 8: 3D Plot of the relative strength i? ^ of Type II versus Type I seesaw mechanisms 
on the complex x plane showing the zeros od R in the region |a;| > 3. 


we are generous with the coefficients (~ 3) that occur in the relation, a 2 cannot decrease 
below 10“^ : so we must maintain at least such distance from the zeros of a 2 so that: 
\a 2 \ > 10“^. The remaining zeros of R are just those of g 2 - We shall see below that 
approach to these zeros is also limited by the USMPs. Thus R cannot viably decrease 
below Rmin = Min{R{xi),i = the values of R at closest permissible approach to its 

zeros Xi- 

No further survey of the parameter space is thus required than to evaluate Rmin- 
To estimate Rminixf) is straight forward given our established USMP-SMP numerical 
techniques pm [T]. We allow an approach to the zeros of R as closely as will permit 
0(2 > 10“^. This limits the radius of the circle of closest approach to the zeros x = 
{1/3,1, (0.198437, —.0992186 ± 2.24266i}} of 02 ,for the 3 real zeros and the pair of com¬ 
plex zeros respectively, to be e* = 10“^{.7, 3.9,4., .53}“^ = (.014, .0025, .0025, .02}. We 
observe values of 


{Xiy 6j, Rmin )) {USJVlPs^ } 


(1, .014, .007, (Ac ~ -14, Aw ~ -.11}} 

(.3333, .0025, .015, {USMPs OK}} 

(.198437, .0025, 6 x 10■^ [Ax < -2.9}} (33) 

(.198437, .025, .016, (Ax < -2.4}} 

(-.099219 + 2.2426h .02, 2 x 10■^ 


19 












{USMPs OK,Fii ~ 2.4}} 


The closest one comes to the required values of i? ~ 10“^ or smaller is in the case of the 
complex zeros of q'^. However even in this case the value of Ffj is simply too small and yields 
a maximal neutrino mass of about 4 x 10“®eH(since / ~ 10“^). To illustrate the behaviour 
in this most favourable but still unviable case we exhibit the variation of R, Ax, Fjj on a 
circle of radius .02 around a complex zero of gs as Figs. 9,10,11. 

{{X0,Ab3 (k-kO), A,f7,y,Ybar}, {-0.0992186 + 2.24266 1, 0.02, 1., 0.95, 0.7, 0.6}} 



Figure 9: Plot of i? on a circle of radius .02 around Xq = —.099219 + 2.2426f 


Finally we turn to the zeros of q 2 - They are not zeros of a 2 so there is no constraint on 
approach to the zero from perturbativity but the conclusions are equally negative : 


{x, 6, •) {USMFs}} = 

{(3 + iV7)/8, .06, .015, {A^ < -.9, Fji ~ 80}} 

{(3 + iV7)/8, .005, .0014, {Ax < -1.3, Fn ~ 2100}} (34) 

The rapid fall of Ax (which also promotes Fjj growth) as one approaches the zeros 
excludes them long before the required small values of R are reached. The strong decrease 
of Ax for values of x near a; = 0 is obvious from Fig. 2 and Fig. 5. 

To sum up : the zeros of R are ringed by zones of exclusion that ensure the 
parameters of the MSGUT can never be chosen to ensure Type II domination 
while maintaining viable USMPs. 
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3.2 Type I fails : Collapse on the final lap. 

We next turn to a consideration of the overall magnitude of the Type I seesaw masses which 
we have shown to dominate the Type II contribution almost everywhere. The Type I mass 
formula reads -up to an irrelevant overall phase- 


M/ = 1.7 X IQ-^eVFi h ^ j3 

p _ 19 ((l-Sx) / 

^ 2y/2r]\\ x{x‘^ + 1 ) ^ 

= (35) 

The complicated effect of the RG thresholds and evolution is all contained inside A^. 
However this quantity is itself limited to lie above —1 by Baryon lifetime > 10^^ yrs! 
Before actually inspecting the effect of Ax in detail we should therefore inspect whether 
the almost completely known function Fj could ever achieve the magnitude of 10 or so it 
would require to yield neutrino masses of around .05 eV (recall h ~ .3 and Ax > —1) in 
the most favourable case. This seems a rather mild requirement but we shall show that 
it cannot be satisfied viably anywhere in the parameter space of the MSGUT 
! 

The limiting behaviour of |R/|is 


x| oo ^ \Fi\ —.02 

x\ —>• {0, ±i} ^ \Fi\ —oo 


(36) 


The approach to the asymptotic value is rather rapid while the singularities as a; ~ 
{0, ±i} are quite mild : ~ |Aa ;|“2 . Figs 12,13 which are 3D plots of \Fi\ over the complex 
X plane for |x| < 10 (more than sufficient to reach the limiting form) and |a;| < 1.5 (to 
show the singularities at a: = 0 , more clearly). 

The dependence of Fj on the ‘slow AM’ parameters { 7 , 7 , 7 , A} is rather marginal- as 
long as we do not consider pathological cases where 7 or A are very small ( forbidden by the 
USMP variation PJ) or pathological cases where perturbativity in any of of these couplings 
is lost. We have verified this by examining 3-dimensional plots of Fj over the complex x 
plane for random complex values of these couplings without hnding any visible variation. 
Apart from the constraint that there magnitudes are greater than .2 their values can be 
chosen randomly without any appreciable effect. 

In Fig. 12 and Fig. 13 the rapid approach to the asymptotic value of |F/| beyond 
|x| = 1 . 2 , the ridge associated with x = 0 (due to its indeterminate phase), and the 
pointed ridges associated with x = ±i, as well as the boundedness by about .1 except 
right close to the singularities, are the relevant and crucial features. Plots which have very 
different coupling values are almost indistinguishable from Figs. 12,13. This is easy to 
understand in terms of the expression for Ziq (Table II) since that is where the dependence 
on these parameters enters. It is evident that the variation of the slow parameters has 
no appreciable effect at all. So they can be set to have any values with magnitudes ~ 1 
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[X = l , r ] = l , ybar = l} 



10 


Figure 12: 3 D Plot of \Fj\ on the x-plane. 


without any loss of generality. The behaviour of F implies that only very close to the 
singular points x = { 0 , ±i} is there any hope that |F/| will achieve the required values 
> 10. Elsewhere it is manifestly bounded by .1 or so and is thus a factor of 10^ short of 
the required magnitude(in the most favourable case of Ax=- 1 : which value is itself only 
borderline compatible with Super-Kamiokande data on baryon lifetime ) 

In Fig. 14,15 we exhibit contour plots of Fj and F/respectively on the complex x-plane. 
Evidently only around x = 0,±i is there any chance of large enough Fj However neither 
region gives a viable result. Firstly, it is evident from Fig.2 that the region around a; = 0 
is forbidden since Ax < — 1. A ‘scissors’ operates between the exploding effect of Ax as 
one approaches x = 0 and the rapidly decreasing \Fi\ as one recedes from it. 

Around x = ±i, Ax is > 1 and hence will suppress rather than enhance the favourable 
value of F/. This is evident in Fig. 15 where x = ±i are out of the running altogether. 
One is again caught in a scissors . At this point (of enhanced G 3211 symmetry) Ax is large 
and positive at the singularity. So it completely overwhelms the growth of Fj due to the 
singularity. 

Thus we have shown that even via Type I seesaw there is no possibility of 
achieving neutrino masses larger than about 5 x lO^^eH in the MSGUT !. 

Although the argument given above is complete we also studied the behaviour at other 
special points such as a; = { — 1,1/3, .5} and the zeros of {q 2 {x), q'^,psi^x),p^i^x)}, some 
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of which have been claimed (and some not !) to represent candidate points for Type I 
enhancement j4dj but are not in fact so. We did not hnd any appreciable growth of Fi with 
tolerable values of Ax at any such point : exactly as expected from our arguments given 
above. Some growth is observed for the x values (digits truncated ) 0.1237 (real zero of 
Ps), 0.33333 (enhanced symmetry point with symmetry SU{5)fupped x 17(1)) and .198 (real 
zero of ^3 ). In all three cases it is due to catastrophic fall in Ax. The contour plot of 
Fj neatly summarizes all this behaviour since only the forbidden region around x = 0 ever 
gives large enough F/. 

To trace the reason for the divergence in our candidates from those ofjlHl is easy : if 
we translate the expressions of |12] which are written in terms of m\m^ — m} into our 
expressions written in terms of h, / we immediately obtain that the functions //, fu of m 
need to be multiplied by (in our notation) {r3/(4r2),r 2 } rsj {(PsPi)/{P 5 P 2 ), {P 5 P 2 )/(psPi)} 
respectively and this then removes the aforementioned points from the list of candidates , 
leaving , after all repeated factors are cancelled, only {0, ±i} as candidates for strong Type 
I seesaw and the zeros of R namely x = {1/3,1} and the zeros of q 2 and q'^. 

4 The Next MSGUT 


We have seen that the Type II seesaw mechanism does not work when only the 10 and 126 
FM Higgs irreps are employed and even the Type I mechanism fails even though not as badly 
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ContourF|Ploli(Cont->.01, ?j= 0.9^ y= 0.7jgambahd0.6> 
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-3 -1.5 0 1.5 3 

Figure 14: Contour Plot of \Fj\ on the a:-plane with contours at \Fj\ = .1, .9. Lower values 
are shaded darker. The small white islands are around a: = 0, ±L 


as Type II. It is natural then to look for a role for the remaining FM Higgs multiplet type 
allowed by SO(IO): the long neglected 120-plet irrep. This irrep has a Yukawa coupling 
matrix that obeys qab = —Qba he. is family index antisymmetric. As such it introduces 
novel properties and problems into the fitting matrix problem. It has so far been cast in a 
minor supporting roleJElEl] in the story of the seesaw. Only a few ant hors im UHl CHI CH] 
have considered the 120plet in a large role and even they did not analyse the htting 
problem comprehensively. The 10 and 120 irreps may together successfully explain the 
charged fermion spectrum, particularly since the 120-plet also contains,besides a (1,2,2) 
sub-representation, a (15,2,2) subrep (w.r.t the Pati-Salam subgroup). Thus it can also 
be used to implement -albeit in a novel way due to its antisymmetric Yukawa couplings- 
a Georgi-JarlskogjHH] mechanism to explain the approximate equality of mg and at 
Ms in the MSSM. This relieves the 126 of the multiple loads it has shouldered so far 
and which, hnally, the present papers show it tends to cast down on the last lap of tests. 
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Figure 15 : Contour Plot of \Fj\ on the x-plane with contour at \Fj\ = 10 lower values are 
shaded darker. Note that the required large values can only lie right close to the forbidden 
region around a; = 0 where Ax < — 1. 


Once the 126-plet coupling is relieved of the necessity of being comparable to the 
couplings of the 10-plet the Type I seesaw mechanism can receive the required further 
10-100 fold boost in magnitude since the Type I seesaw masses are proportional to f~^. 
It should be emphasized that the smallness of / in this scenario must be considered as a 
structurally dehning condition rather than as the smallness of a perturbation because it 
leads to light right handed neutrinos. A similar scenario was considered in^H] but analyzed 
by making various somewhat arbitrary assumptions. Thus no generic analysis like that for 
the 10 — 126(21112312211211 system is vet [IH] available to compare with. 

In this way we propose that our results provide a hint of the direction to proceed: every 
member of the SO(IO) FM Higgs family will receive its characterstic function, role and 
gainful employment. It seems that the Good God may be just at least to the different 
breeds of Higgs particles! 
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To implement our (long feltjlH]) intution of a proper share and role for each Higgs type in 
the intricately wrought elegance of SO(IO) unihcation we have calculated jlHl the complete 
set of couplings and masses of the Next to Minimal Supersymmetric GUT or NMSGUT 
defined by the inclusion of a 120-plet Higgs into the MSGUT Higgs complement. Here we 
restrict ourselves to a brief survey of the essential features. 

In the notation of|2D] the real vector indices of the upper left block embedding of SO(6) 
in SO(IO) are denoted by a, 6 = 1,2..6 and of the lower right block embedding of SO(4) 
in SO(IO) by d,/3 = 7, 8, 9,10. The doublet indices of S'f/(2)L(S'f/(2)ij)are denoted by 
a,/9 = l,2(d/9 = 1,2). The index of the fundamental 4-plet of SU(4) is denoted by a 
(lower) yU, i/=l,2,3,4 and its upper-left block SU(3) subgroup indices are /2, h = 1, 2, 3. 

The decomposition of the 120-plet w.r.t SU{A) x SU{2)l x SU{2)r is as follows. 


Oi,fc(120) = OaU^O + TO, 1,1) + Oabd^b, 2,2) + 1,3) + (6,3,1)) + 0^^.^(1,2,2) 

= OW(10,1,1) + Of;,(IO, 1,1) + 0„„i'‘(15,2,2) 

+ 0'“j.y6,l,3) + 0<“yy6,3,l) + 0„i(l,2,2) (37) 

(where we have used the superscripts ^ discriminate the symmetric 10-plet from the 
antisymmetric 6-plet). 

The decomposition of its coupling to 16-plets of SO(IO) is jSOl 

+ (38) 


The SU{4) X SU{2)l x SU{2)r labels in the above jHII] equations are trivially decomposable 
to SM labels and thus such expressions - available for every coupling in the MSGUT and 
NMSGUT-constitute a complete solution for the Glebsches to these MSGUTs. 

The 120-plet contributes two pairs = 0“, = O^, hg = 

of MSSM ([1,2,1]0 [1,2,-1])doublets. Their Yukawa couplings to MSSM fermions follow 
immediately from the above unitary decompositions: 


= v{h + f + g) ; ri = — cot (3 ; r 2 = — cot f3 

CKi 0(2 




= v{{h - 3/) + (rg - 3)^) ; rg = 


m = 


m = 


4i\/3a5 

n(rih-l-r 2 / + r6^); rg = cot (3 

CKg T 1y dCKg 

n(rih-3r2/-f-r7^); fs = cot [3 

ttg + '^y dctg 
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2igJ-{ae + iVSa^) sin (3 ; r 7 = (rg - Sre) 


(39) 


which clearly exhibits the analogy with the Georgi-Jarlskog type couplings in the MSGUT. 
We propose to fit the charged fermion masses using light doublet components of predom¬ 
inantly 10 and 120-plet origin. The additional terms in the superpotential do not mod¬ 
ify the AM SSB and moreover still contain only the 26 MSSM multiplet types that we 
describedjSI] for the MSGUT. The corresponding mass matrices- some of whose dimen¬ 
sions, like the doublet mass matrix given below, are raised relative to the MSGUT- will 
be given, in our notation, and along with the corresponding RG and USMP-SMP analysis 
inj46j. A corresponding cacluation of mass matrices only - with very different conventions 
- is already available jUj. 

The extra terms contributed by the 120-plet to the superpotential are : 


lUl20 


2^21^ T Oijk^ijlmn^klmn T '^J^^^^Oijk^ijlmn^klmn 


(40) 


The coupling of the P([3, 3, ±|])and K([3,1, ±|])-type color triplets contained in the 
120-plet leads to quite novel, family index antisymmetric and even SU{2)^ triplet mediated, 
contributions to the effective LLLL and RRRR type AB ^ 0 violating d=5 operators in 
the low energy effective theory |46j. These may lend additional texture to the emerging 
story of the deep connections between neutrino mass and Baryon violation j49[ l!?T] . 

The composition of the MSSM [1, 2, ±1] doublet sector can be deduced from the left 
and right null eigenvectors of the 6x6 mass matrix of the NMSGUT analogously to those 
dehned for the MSGUT in Section II. 


/ 

-Mh 

'Y\/3(uj — a) 

-7\/3(i.i4 + a) 

— 7cr 

kp 

— y/SikiO 

\ 


-7-\/3(aj + a) 

0 

— (2M + 4r](a + cj)) 

0 


i{p + 2ij)C 



^y/Z{LO — a) 

— (2M + 4ri(a — uj)) 

0 

— 2r/d-\/3 


-i{p- 2uj)C 



—a'j 

—2r}u\/li 

0 

— 2m + 6A(Lt; — a) 

Co- 




pk 




-Mo 

-^iuj 

\/3 


V 

\/3ikLj 

3 

1 

-i(p + 2uj}( 


— 

-Mo - 

/ 


The left and right null eigenvectors are calculated after imposing the light doublet 
condition DetTi = 0. 

The relative strength of the contributions, of the 10 and 120-plets vis a vis the 126- 
plet, to the charged fermion mass matrices can be studied by examining the variation of 
« 2 /tti, 5 , 6 , d; 2 /d;i, 5,6 in exactly the same manner as already done by us in this paper. The pure 
form of this scenario can thus be examined directly in terms of htting the charged fermion 
masses and mixings in terms of 10 and 120-plet contributions only. If this succeeds to a 
within about 10%, the effect of including the 126-plet contributions can be examined j46j. 
Naturally, the investigation is best begun with the 2x2 case as for the 10 -|- 126 system 
pn |2ni CS] • This is reported in |17j . 








5 Conclusions and Outlook 


The Supersymmetric GUT based on the the 210 © 126 © 126 AM Higgs system is the 
simplest Supersymmetric GUT that elegantly (almost, v.supra!) realizes the classic pro¬ 
gram of Grand Unification. Its symmetry breaking structure is so simple as to permit an 
explicit analysis of its mass spectrum at the GUT scale and an evaluation therefrom of the 
threshold corrections and mixing matrices relevant to various important quantities. It im¬ 
plements in a fascinating and elegant way an intimate (“ouroborotic”) connection between 
the physics of Lepton number and Baryon number violation. Since it has the least number 
of parameters of any theory that accomplishes as much this theory merited the name of the 
minimal supersymmetric GUT or MSGUT. The same simplicity and analyzability of GUT 
scale structure also applies to the theory with an additional 120-plet, since it contains no 
standard model singlets, and thus justifies calling it the Next to Minimal (or New Minimal) 
Susy GUT (NMSGUT). The small number of Yukawa couplings of the MSGUT makes the 
fit to the now well characterized fermion mass spectra very tight. Inin we initially ob¬ 
served and in this paper we have shown that the combined constraints of the 
seesaw fit and the preservation of MSSM one loop gauge unification are enough 
to rule out the MSGUT. 

We now argue that the very natural inclusion of the third possible SO(10)FM Higgs 
type, i.e.the 120-plet, which was somewhat arbitrarily [H] -but very fruitfullyUni EHl 113 El 
ini El 122 El El 123 -excluded from a leading role in the ‘SO(10)party on the GUT scene’, 
offers a natural resolution of the difficulty of the MSGUT in achieving large enough Type I 
seesaw masses. In this scenario the 10-plet and 120-plet are primarily responsible for the 
charged fermion mass £t due to a relatively suppressed contribution of the 126-plet to the 
charged fermion masses due to a small value of the Yukawa coupling 16 ■ 16 ■ 126. This very 
suppression implies that the Type I seesaw will become enhanced essentially due to lowering 
of the righthanded neutrino masses. Thus far from being a perturbation the role of the 126 
is rather to define a very specific class of Grand Desert in which the right handed neutrinos 
have masses considerably less than the generically expected GUT scale masses. Thus this 
proposal, which was thrust upon us by detailed analysis of compatibility between GUT 
and Neutrino mass scale hierarchies, is distinct from previous scenarios and, moreover, has 
very clear and distinct phenomenological consequences. Furthermore the 126-plet Yukawa 
coupling is itself released from the Type I charged fermion constraints and is subject only 
to the relatively mild constraints coming from the right handed Majorana neutrino mass 
limits from cosmology. If this does not conflict with the requirement to achieve the small- 
large /quark-lepton mixing duality then the one degree of magnitude failure of the seesaw 
in the MSGUT should be quite surmountable. In a slogan: 

The MSGUT is dead! Long live the (Next) MS GUT! 
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7 Appendix 


In this Appendix we express our formulae [311 |T] for neutrino masses etc in terms of the 
variable x which parameterizes ma the fast variation of AM scale mass matrices in an 
elegant and simplifying manner, as has been particularly emphasized bv|321 1^ . This 
allows the easy analytic localization of putative points of seesaw growth|43j. It enhances 
conhdence in the completeness of the exclusion by survey PP of any possible seesaw £t. The 
use of the variable x also signihcantly simplihes the presentation of data since considering 
variation over the x-plane unihes consideration of the three solutions for x corresponding to 
a given value of ^ while retaining complete clarity due to the rational (3 to one) map from 
the x-plane to the ^ plane. The process of translation merely makes explicit in the mass 
matrices what was already done[TniinillS21 when parameterizing the AM vevs in terms of x. 
If we substitute for the AM vevs in terms of x, the equations for the Null eigen-vectors oi Ti¬ 
the doublet mass matrix of the MSGUT - subject to the condition DetTi = 0 become: (note 

the labelling of doublets h®h= (i/f , ^)) 


A 

A 

A 

A 

N 

N 


0^2: *^3? CX 4 } — ^{^1') ^2') *3^3, 0 ( 4 } — NA 
0^2: ^3: *^ 4 } “ ^2? *^3? ^ 4 } “ NA 

(737^4(1 -X)) (737^2(1-2:)) 7 g 3 d(x-l) 

2r]p5 ’ 27 P 3 ’ 2a/^P5 

(\/37P4 (1-a;)) (\/3 7P2 (1 - a;)) 

2r]p5 ’ 2pp3 ’ 2^/rjJp5 

\/©6 

IPsPbl 


(41) 


The crucial functions Fj, Fjj, R used as SMPs[T] in our survey now become (up to 
irrelevant overall parameter phases ) 
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Fi 

Fn 

R 


10 7C/ I I (1 - 3x) q'^ 

2V2 + 


^q-Ax 2 ^ 27 ^ \p 2 P^Pb\ ^ + 1) (4a; - l)g| 

y/^ (x - 1) V zie \ x(l - 3x) q3q2P5 

I Fi ^ (x- l)(3x- l)g293^ | 

Fjj 8(4x — l)(x^ + 1)^3 


(42) 


The various polynomials defined in 1221113 together with some additional definitions 
are given in Table II. 

A number of remarks are in order 

•The asymptotic behaviour of {Fj,Fn,R } as |x| —> 00 is {|x|°, |x|“^, |x|} and as 
X —0 is {|x|“ 2 , kl°} so it is obvious that only a bounded range of (say) |x| < 10 

need be considered to check if R is small or Fj is large. 


• For Type II not to be dominated by the Type I implies the ratio R should be very 
small (~ 10“^). So it would seem that the enhanced symmetry points x = 0, 1/3 and the 
additional points x = {0.198437, -0.0992186 ± 2.24266 i}, which are the zeros of 

g 2 , may help to strengthen the Type II seesaw and should be added to our list of excep¬ 
tional points. If one extracts a factor to define R' = Ra^^ one sees that only the the 
zeros of q 2 are of any relevance to stengthening Type II seesaw . Our different parametriza- 
tion of the problem of dominance together with the extraction of the hidden dependence 
on the threshold corrections to Mx which is coded in Ax makes our analysis jl] different 
from|121- Differences fromj3H| in the points relevant for Type II have emerged: we exclude 
X = 0 and the zeros of P 5 ,P 2 since they cancel with the normalizations (which we have not 
left implicit as in[22])- Moreover we have used expressions in terms of the actual Yukawa 
couplings rather than mass eigenvalues which can conceal compensating growth of Yukawas 
and coefficients necessitating separate consideration of whether one has ventured into 
strong Yukawa coupling regions. Nevertheless we have included all such putative special 
points points in our survey and checked that they do not provide the required exceptional 
behaviour. 

•The functions 2 : 2 , have no zeros even when 7,7 are both zero. When 7,7 are non 
zero the Normalizing factors never diverge. 

•The points x = 0, are the only remaining candidates to strengthen Type I and are 
already included in the list of enhanced symmetry points. 

Note Added : Our surmise was announced in May/June 2005 (PLANCK05,Trieste, 
May 2005 and hep-ph/0506291) and our proof appeared as hep-ph/0512224vl). In May 
2006 a calculation appeared PI confirming that optimal FM fits found by the “downhill 
simplex” method are not viable in the MSGUT, for exactly the reasons given by us. 
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— 3a; -I- 1 


4a;^ — 9a;^ -|- 9a; — 2 


P 4 /(3a; — 1) = a;^ -1- 5a; — 1 

P2 

( 2 a;-l)(x + l) 

P3 

12a;=^ - 17a;^ -flOa; - 1 

Pi 

(3a; — l)(a;^ -f- 5a; — 1) 

Pb 

9a;^ -I- 20a;"‘ — 32a;^ -|- 21a;^ — 7x -|- 1 

Z2 

2 a; p -f- 1 — a; p 

Zie 

In.-nuP 1 3|7G-1)P3P4|2 , 3 1 7G-1)P2P5 |2 , 1 \2 

mPbl 1 4 I r, 1 ' 4 I v 1 ' 2 !/'^^ ^jPsQsy Xr, \ 


\r>.r)W 1 3 | 7 G- 1 )P 3 P 4|2 , 3 1 7 G- 1 )P 2 P 5 |2 , 1 1^7 .x Jx{l-bx){l+x'^) 

^16 

\PWb\ 1 4 I ^ 1 1 4 I r) 1 ' 2\’W ^jPsqsy Xr, \ 

N 

IPsPbl/x/^ 

N 

IPsPbl/V^ 


Table 2: Functions of x entering the expressions for the USMPs and SMPs. 
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